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We propose comparing cosmological solutions in terms of their total spatial volumes V(t) 
ON , as functions of proper time r, assuming synchronous gauge, and with this intention evaluate 

ON ■ the variations of V(r) about the Friedmann-Lemaitre-Robertson- Walker (FLRW) solutions 

for dust. This can be done successfully in a simple manner without solving perturbation 
\ equations. In particular, we find that first variations vanish with respect to all directions 

which do not possess homogeneity and isotropy preserving components; in other words, every 
FLRW solution is a critical point for V(r) in the properly restricted subspace of the space of 
solutions. This property may support a validity of the interpretation of the FLRW solutions 
as constituting an averaged model. We also briefly investigate the second variations of V(r). 



' 81. Introduction 

co ; 

C The most important practical models in relativistic cosmology are the so-called 

Friedmann-Lemaitre-Robertson- Walker (FLRW) models, which assume the homo- 
geneity and isotropy of spatial geometry and matter distribution. Since our universe 

ON , seems homogeneous and isotropic over a large scale (typically the size of a superclus- 

ter), it may be natural to expect that such a model provides a good approximation 
for the real universe, describing averaged properties of the geometry and matter dis- 
i ', tribution in it. However, t!P justification of such an expectation should be considered 

carefully, since Einstein's equation is highly nonlinear, and thus no simple averaging 
procedure applied to spatial geometry and matter commutes with time evolution. 

. i-h [ The averaging problem in general concerns how one can modify Einstein's equa- 



tion for averaged variables to be compatible with the time evolution of the original 
inhomogeneous data for the variables. To analyse the problem one usually has to 
determine^naveraging procedure. In fact, many averaging procedures have been 
proposed, B ~i3 but the dynamical properties of the resulting averaged variables are 
often unclear, and little firm agreement is seen among them. 

In this paper we discuss the problem in such an indirect way that no explicit 
averaging procedure is introduced, so that the results do not depend on the aver- 
aging procedure. We simply compare the dynamical properties of the solutions - 
through the temporal behavior of their total spatial volume V(r), and we find a 
striking property of the FLRW solution for dust matter. One of the solutions that 
is compared is an FLRW solution gj^, and the other is an inhomogeneous solution 
gal,. Synchronous and comoving conditions are assumed. The reason we focus on 
the volume V(r) is that the FLRW model has only one dynamical variable, the scale 
factor a(r), and its dynamical content is equivalent to the volume Vq(t). Thus, it is 
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reasonable to extract V(r) also for an inhomogeneous solution and use it to compare 
the solutions. To be specific, one could suppose the norm of the difference between 
the solutions defined by, e.g., the L p -norm applied to the volumes, 

l|g o6 -g2 ) ll P = (/_jnr)-%(r)rdr) 1/P , (M) 

though we will not use any specific norm. The difference AV(t) = V(r) — Vq(t) 
itself is the object on which we focus. 

Our main claim is that each element of the two-parameter set of the FLRW 
solutions for dust is a critical point for V(r), i.e., the variations of V(t) in the 
space of solutions about the FLRW solution vanish: 5V(r) = 0. (This statement 
shall be made more precise.) Due to this property the temporal behavior of an 
inhomogeneous solution which is almost homogeneous and isotropic is the same as 
that of an FLRW solution up to first order in the sense that AV{t) ~ 6V(t) = 
0. This implies that there exists a natural correspondence between inhomogeneous 
solutions and the FLRW solution. We interpret this correspondence as an averaging 
(up to first order). On the other hand the second variations 5 2 V(t) do not vanish 
in general. We will also give some formulae for the second variations. 

The content of this paper is basically a direct generalization of a previous paper, 
where we examined the spherically symmetric case using the exact Lemaitre- 
Tolman-Bondi solution. The method used to prove the claim in this paper is different 
from that in the previous paper, since we cannot rely upon an exact solution. The 
existence of a conserved quantity (the total number of dust particles) will be one 
of the key points of our analysis. Another key is the property of a homogeneous 
and isotropic (three-)metric, mainly as an Einstein metric, i.e., that the Ricci tensor 
is proportional to the metric. These two features give rise to a closed ordinary 
differential equation for 6V(t), and which is used to show the claim. In particular, 
we do not have to solve perturbation equations. The second variation 5 2 V(t) is 
obtained in a similar way. 

Our approach is not to present a complete framework for averaging, but rather 
to give a transparent view of the interrelation between inhomogeneous solutions and 
the FLRW solutions. This should be a helpful guide for a detailed investigation. 

In the next section we make some preparations for our analysis and introduce 
some notation. In §|3| we discuss the first variations of V(r) mainly in the case of 
closed spatial manifolds. Section |I| concerns the second variations. In §||we comment 
on cases of spatial manifolds which are compact but have boundaries. Section || is 
devoted to a summary. 

§2. Preliminary 

Let M be a 3-dimensional compact manifold which admits a constant curvature 
metric hij. We use x l to denote coordinates for M. The spacetime manifold we 
consider is the direct product M x R, for which we use x° = r to denote the 
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coordinate for R. We consider smooth metrics on M x R of the synchronous form 

ds 2 = g ab dx a dx b = dr 2 - hijir, a^dzW. (2-1) 

Here the spatial metric hij(r,x l ) can be considered as a one-parameter (r) family 
of smooth metrics on M. Conversely, once given a one-parameter family hij(r,x l ) 
of smooth metrics on M, we can identify it with the spacetime metric fl24|) . Thus, 
we can define the space P* of smooth (synchronous) metrics on M x R as the set 
of all possible one-parameter families of smooth metrics on M, /iy(r,ic l ). Next, let 
us consider Einstein's equation for dust: G a b = 8-KpUaUf,, where G a b is the Einstein 
tensor for g ab , and p and u a = (d/dr) a are, respectively, the matter density and 
four- velocity of dust. We define the space of solutions P C P* as the set of all 
smooth solutions of this equation. We distinguish P and P* to clarify whether a 
relation holds with dynamical equations or without them. 

Since we have assumed that M admits a constant curvature metric hij, the 
spacetime manifold admits a spatially homogeneous and isotropic metric of the form 

ds 2 = gf b ] dx a dx b = dr 2 - a(r) 2 %(s ! )dxM^. (2-2) 

The function a(r) specifies the scale of M for a slice r =constant, and it is called the 
scale factor. To make the magnitude of a well-defined, we remove the scale ambiguity 
of hij by choosing it (IS i\- standard metric, satisfying 

(3) R = 6k, (2-3) 

where is the scalar curvature for the standard metric chosen, and k = 1 for 
a positive curvature space, k = —1 for a negative curvature space, and k = for 
a flat space. The constant k is referred to as the curvature index. The curvature 
tensor ^Rijkh Ricci tensor ^Rij, and scalar curvature ^R for the spatial metric 
hij = a 2 hij can be found, for example, from the conformal invariance of the Ricci 
tensor ^Rij'. 

^Rij = ^Rij = 2khij = -phij. (2-4) 

The curvature tensor and scalar curvature are hence 

k 6k 
{3) Rijki = -^{hikhji - h u h jk ), ®R = -p. (2-5) 



We define the total spatial volume V(r) for the metric (24) by 



V(t) = / \ r hd 3 x, (2-6) 

JM 

and the total particle number E by 

E = f G 00 Vhd 3 x = - f ( ( % + K 2 - KalC^y/h^x, (2-7) 
Jm 2 Jm 

where ^'R is the scalar curvature for the spatial metric hij, Kij = \hij is the extrinsic 

curvature, and K = h l] Kij, h = det(/ijj). To raise and lower the spatial indices 



4 



M. Tanimoto 



i,j,--- we use and hij, e.g., K l j = h tk Kkj. Here dots (') represent derivatives 
with respect to r. 

The particle number E is conserved if the dynamical and constraint equations 
Gij = = Goi are imposed. In fact, one can check by a straightforward calculation 
that (G oVhJ = (-GijK ij + V^o^Vh, implying E = if we impose the equations. 

V(t) can be regarded as a "function- valued" functional on P (or P*), since once 
hij(r,x l ) £ P (or P*) is specified the function V(r) is determined. (To express this 
we could add an argument of the functional as "V(r)[/iy(r)]," but for simplicity we 
shall write "V(t).") Similarly, P is also a functional on P or P*, which is real- valued 
on P, since E 1 is conserved, but is function- valued on P*. 

Consider a smooth path I in P*, I : [0, oo) — > P*. Let K(r) be the function 
V(t) for i(e), where e G [0, oo). Then, we can expand 14 (r) at e = as 



K(r) = V (r) + e<5V(r) + ^e 2 8 2 V(r) + • • • , (2- 



where 



de 



(2-9) 

=o 



We call 8V{t) and 5 2 V{t) the variation and second variation of V(r), respectively. 
Variations of any functional in any space are defined similarly. 

Practically we will not specify any path, but we expand a variation of a functional 
in terms of variations of hij and Kij. Thus it is useful to introduce some notation 
for them. 

Notation Let Shij and 6 Kij be variations of the spatial metric and extrinsic cur- 
vature for the spacetime metric (2-1). Similarly, let 5 2 hij and 5 2 K{j be their second 



variations. We use the following notation for them and their traces throughout this 
paper: 

lij = Shij, 7 = Xij = 5K ij7 X = h lj Xij, (2T0) 

j^ij = 6 2 hj, 7 ( 2 ) ^ W^ij, A^,, ^ 5 2 K tJ , \W = h ij X ( 2 hj. (2-11) 

§3. Variations of V 

Our concern is the variations of V(r) and E evaluated about the homogeneous 
and isotropic metric ( |2-2j) . We find a closed relation between them if the spatial 
manifold M is closed, but for later convenience we present the relation in the case 
that M is compact first. 

Lemma 1 Let V{t) and E be the junctionals on P* defined by Eqs. (\) and (p?- 



For the variations about a spatially homogeneous and isotropic metric ggh , 




where 



SE - aSD j — .", ( - ) )5V + 2^5V, (3-1) 



8B = ±- [ V^V^ij - V a )Vhd 3 x, (3-2) 
2o Jm 
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a = a(r) is the scale factor for the homogeneous and isotropic spatial metric, and k 
is the curvature index. 

Here, 5V is the time derivative of 5V ("5" and " ' " commute, however), and Vj 
represents the covariant derivative associated with (the zeroth order of) hij. As long 

as there is no confusion, we omit superscripts (such as that on hf^) in reference to 
a zeroth order quantity. 

To prove Lemma |l| we present some formulae, which can all be checked by 
straightforward calculations. 

Formulae 2 The variations 5V , SV , and 5E (about a generic point) are given by 



SV = - 



2 Jm 



1 



Vhd 3 x, 



SV = j (A + K<y - Kijf^ Vhd 3 



(3-3) 
(3-4) 



and 



5E = J - Q + KK ij - K\Ki k 



+ 



+ K 2 - KijKV) 7 - (K ij - Kh ij ) A 



+ 2 Vi ( V ^'-V 47 ) 



Vhd 3 x. 



(3-5) 



Proof of Lemma [J. We evaluate the variations fl3-3j ), fl3-4|) and (3-5) about a 
homogeneous and isotropic metric. The spatial part of such a metric can be written 



(3-6) 



as hij = a (r)hij, as in Eq.(|2-2|). The extrinsic curvature becomes 



Ki 



2 hij 



a 



The Ricci tensor ^Rij is also proportional to hij (Eq. (|2-4f )). If we substitute Eqs. (|3-6|) 
and (p]l) into Eqs.(|l|) and (3-5), we find 



and 



6V 



5E - a5B 



M 



la 



A + --7 I Vhd 

2 a 



x 



X^hd 3 x + -5V 
M a 



) / jVhd 3 x + 2- / Av^d 3 ^ 
J Jm a Jm 



(3-7) 



(3-8) 



The last equation (3^8 ) implies Eq.(3d) if we note 
0), the divergence term 5B 



where SB is defined by Eq 
Eqs.Q and Q. □ 

If the spatial manifold M is closed (i.e., dM 
vanishes, so that we have a simpler result: 

Corollary 3 Let V(r) and E be the functionals on P* defined by Eqs.(2-t) and 
( [j?- / | ). // the spatial manifold M is closed, for variations about a spatially homoge- 
neous and isotropic metric g^, 



5E 



SV + 2-SV, 
a 



(3-9) 
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where a = a(r) is the scale factor for the homogeneous and isotropic spatial metric, 
and k is the curvature index. 

The relation (3^9) also holds for the variations taken in P. In addition, we can 
regard 5E as a constant, since E is conserved in such a case. The scale factor a(r) 
is also a definite function which obeys the Einstein equation 

d ( d\ 2 k , 
2-+ - + ^ = 0. 3-10 
a \aj a z 

Hence the relation ( |3-9|) is a closed linear ordinary differential equation for SV. 
Although we can easily find the solution for Eq. Q3-9| ) directly, for later convenience 
we present a more general formula first: 

Formula 4 The general solution for the ordinary differential equation for f = f(r) 

g(r) = (^-3(^) 2 )f + 2^f, (3-11) 
\ or \aj a 



where g(r) is a given function of t and a = a(r) is a solution for Eq.(S-lC), is given 
by 

.2,(1 faij). 

,2 7 ad 2 



f(r) = a 2 d{- I ^dr + c), (3-12) 



where c is an integration constant. 

Applying this formula we at once obtain the following result. 

Lemma 5 Let V(r) be the functional in the space of solutions P defined by Eg. ( \2-(\ ). 
If the spatial manifold M is closed, variations about a homogeneous and isotropic 
solution are given by 

5V{t) = a 2 d J ^ + SCj , (3-13) 

where 5E and 5C are constants, and a = air) is a solution of Eg. ( p-TZj j. 
We present the explicit forms of the solution for completeness: 
(i) k = 1: o(r) = a (l - cos??), r - r 0c = a (rj - sin??), 

SV{t) = Oq(1 — cos if) (A{rf)-^- + sin rjSC) , (3-14) 



where A{rf) = (1 — cos??) 2 — sinr?(?? — sin??) = 2(1 — cos?/) — ??sin??. 

(ii) k = —1: a(r) = ao(cosh?? — 1), r — tq c = ao(sinhr? — r?), 

SV(t) = a 2 , (cosh?? - 1) (a.-(v)-^- + smhr]6c\ , (3-15) 

where A-{rf) = (cosh?? — l) 2 — sinh ??(sinh ?? — rj) = r/sinhr? — 2(coshr? — 1). 

(iii) k = 0: a(r) = a (r - r 0c ) 2/3 , 

SV(t) = \{r~ roc) 2 SE + ^a 3 (r - r 0c ) SC. (3-16) 

In these solutions ao and tq c are constant parameters. As one can easily see, t$ c 
is redundant as far as the FLRW solutions are concerned, since it carries the gauge 
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freedom associated with the choice of the origin of time. However, we will see that 
roc is of some importance in the wider context in which we are interested. 

Does the constant 5C in Eq.( |3-1 3] ) defined as an integration constant for the 
differential equation ( |3-9| ) have an "integral" C such that its variation coincides 
with the SC that connects variations of E and V(t) through Eq.( [3-13| )? We may 
expect that at least in a neighborhood of the FLRW solutions there exist such a 
single-valued functional C . We will call this functional C, a big bang constant. For 
example, in the space of spherically symmetric solutions there exits a natural choice 
of C that is globally defined, 03 but in the general case, global existence of C, of 
course, requires proof. However, we do not discuss this problem further, since local 
existence is sufficient for our purposes. We do not discuss the uniqueness of C either, 
for the same reason. (In fact, it would not be unique from the definition above. A 
complete definition may be given after higher order variations of V(r) are specified.) 
We will see that considering the functional C is convenient for our analysis. 

The function V(t) for an FLRW solution, denoted Vq(t) hereafter, depends 
parametrically on ao and tq c , and thus it varies if ao and tq c are varied (with r 
fixed). For example, in the case k = 1 with M ~ S 3 , the total derivative dVo with 
respect to ao and tq c is0) 

dVb(r) = 67r 2 aQ(l — cos ri)(A(rj)dao — sinr/dro c ). (3T7) 

Note that the time dependence of this function is the same as that of Eq. fl3T4|) , 
Similar results hold for the other cases of k. To understand the implication of this 
fact, we first need to give some definitions. 

Definitions Let O be the two-dimensional space of the FLRW solutions, and let 
o(ao,To c ) £ O be the FLRW solution with given FLRW parameters ao and ro c - O 
is a subspace of P, as well. Let T a0tTQc (P) be the tangent space at o(ao,To c ) £ P- 
The tangent space at o(ao,To c ) £ O can be considered as a subspace of T a0)TOc (P), 
and we denote it O a0!Wc (P). We call O a0jTOc (P) the tangent space which preserves 
the homogeneity and isotropy. 

In correspondence to O a0j7J)c (P) we can define a complement Q ao , TO c(P) D Y 
T ao ,TX)c(P) = O ao ,T 0c (P) © Qa ,r 0c (P), where is the direc t sum . Q ao ,T 0c (P) is not 
unique at this point, but a comparison of Eqs.(3T4) and fl3T7| ) reveals that a tan- 



gent space at o(ao,ro c ) spanned by independent tangent vectors with nonvanishing 
5E and 5C is naturally identified with the tangent space O aQ)TQc {P). This identifi- 
cation makes Q ao ,T 0c (P) unique, since we can require that for any -J— £ Qa ,T 0c (P) 
the corresponding variation of V(r) vanish, 5V(t) = 0. To summarize, we have: 
Theorem 6 Let V(r) be the functional in the space P of solutions defined by Eg. \2- 6Q , 
and suppose that the spatial manifold M is closed. Then, the tangent space T aQtTOc (P) 
at o(ao,Toc) £ P can be uniquely decomposed as the direct sum 

Ta ,r 0c (P) = Oa (h r 0c (P) © Qa ,r 0c (P), ^-18) 

where Q ao ,T 0c (P) is the space for which the action of any -J— £ Q ao ,T 0c (P) on V(r) 
vanishes: 

5V(t) = 0, (3-19) 
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and O a0)TOc (P) is the tangent space which preserves the homogeneity and isotropy. 

In the sense of Eq.fljTlSD an FLRW solution is a critical point for V(r) in the 
space of solutions. We may interpret this fact in connection with the averaging as 
follows. Let p £ P be an inhomogeneous solution which is almost homogeneous and 
isotropic, and let V p {t) be its volume as a function of time. Then, the theorem 
implies that there exist the set of values (ao,To c ) for which the difference between 
the volume Vq{t) of the FLRW solution o(ao,To c ) and that of the inhomogeneous 
solution is the same up to first order: AV(t) = V p (t) — Vq(t) ~ 8V{t) = 0. 

This set of values (do, Tq c ) is determined by the conditions 5E = and 5C = 0. In 
other words, the best fit FLRW solution o(ao, to c ) is the solution in O that possesses 
the same particle number E and the same big bang constant C as those of the 
solution p. This correspondence p — * o(ao,TQ C ) defines the map Av : P —* O. That 
is: 

Definitions The space of solutions P is foliated by the level sets for E (i.e., the sets 
E =constant), and it is also foliated by the level sets for C. Let Pe,c C P be the 
double level set E = constant (= E) and C = constant (= C). In Pe,c there exists 
only one FLRW solution o(ao, tq c ). The map Av is defined by this correspondence: 

Av: P E ,c - o(a ,r 0c ). (3-20) 

We call this map the averaging map in the space of solutions. 

Strictly speaking, the map Av can be defined only on the domain where the big 
bang constant C is defined. This domain exists at least in a neighborhood of O. 
Again, this is sufficient for our purposes. 

In terms of the averaging map Av, we can reword Theorem ^ as follows. 
Theorem 7 The variation ofV{r) along any smooth path in Av~ 1 o(ao, ro c ) from 
an FLRW solution o(ao,ro c ) vanishes: 

5V(t) = 0. (3-21) 
§4. Second variations of V 

The variation of V(r) about an FLRW solution is essentially zero in the sense 
of Theorem |6| (or [?]). The second variation 5 2 V(t) therefore gives the leading con- 
tribution to the difference AV(t). In this section we evaluate 5 2 V(t) in a manner 
similar to the previous section. 

We first obtain the following. 



Lemma 8 Let V(r) and E be the Junctionals on P* defined by Eqs. j\2- by and (2-1), 



and suppose that the spatial manifold M is closed. For the second variations about 
a spatially homogeneous and isotropic metric , 

■ \ 2\ 




3(-J I 8 V + 2— <5 V + r(r), (4-1) 



where r(r) is the function of r that is determined from the first variation of the 
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metric, defined by 



r(r) = 



a 7 



M 

+^((Vo)(V t 7)-(V i7 ^)(VV fc ) 



1 



+^((V^)(V J Y fc )-(V j7 )(V i 7 y ) 



+2^ ( 7 * 'Ay - 7 A) + A 2 - 



\fhd 3 x. 



(4-2) 



Here, a = a(r) is the scale factor in gjj9 ' , and k is the curvature index. 

To prove Lemma |8| we need to calculate the second variations 5 2 E and 6 2 V(t) 
and the time derivative 5 2 V, about a generic point in P* or P. We obtain: 
Formulae 9 The second variations 5 2 V, 5 2 V and 5 2 E are given by 



5 2 V = - 



\ J m (l {2) ~ + \l 2 ) Vhd 3 x, (4-3) 



5 2 V 



M 



\W - 2^\ tJ + 7 A - IC* ( 7 ( \- - 27^7/ + 77*; 

Vhd 3 x, 



+\k ( 7 (2) - iHj + \i 2 



(4-4) 



and 



-^ ( %(-7 (2) +7^-^7 2 ) 



4 

+~V i (5,0 < + ^7/5 4 ) 



2 w 2 
+ (K^K kl - K ik K^ 1 ) lijlkl 



(Ki k K jk - KKij) (7 



(2)»j 



2 7 «y fe + 7 yi 
1 

4 V" " "'■'/ V' ' " J ' 2 
(a% - 4 7ifc A/ + 27 -A + 7 A 



+ \ (K 2 - K' J kj j) ( 7 (2) " fHj + J 7 2 



where 



+K (\& - 2^ Xij + 7 A) - A^ Ajj + A 2 



v^d 3 



(4-5) 
(4-6) 
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and 5f3 l is the variation of (3 l = h %3 j3j. 

Proof of Lemma [|. We evaluate Eqs.(4-4) and (4-5) about the FLRW solution, 
using Eqs.(flD, Q, and Q. The result for 5 2 V is 

5 2 V = [ ( A (2) _ 2j ij Xij + 7 A) Vhd 3 x + A / f 7 (2) + 7 ij 7*i - ^ v^d 3 * 

= / f A< 2 > - 2^- + 7 A) v^d 3 x + ^5 2 V 
Jm k J J a 

+ -! (^^-ll 2 )Vhd 3 x, (4-7) 



a Jm V 2 



where in the last equality we have used Eq. fl4-3| ). A straightforward evaluation of 
5 2 E gives 

p E =U±-(i) 2 ) [ 7 ^Vhd 3 x + 2^ [ X^Vhd 3 x 
2 ya z \aj J Jm a Jm 



f 

Jm 



+ 

I M 

+ Iyi( V fc V fc7ij + ViVjj - 2V*V fe7 /) + ^V^V^- - V i7 ) 



-2-7«A<i + (A 2 - 



Vhd 3 x. (4- 



We have dropped the divergence term ~ / Vj(<5/3 4 + ^ 7 /3*)\//id 3 x. Finally, applying 
integration by parts to the terms containing covariant derivatives, and using Eqs.( |4-3| ) 
and (4-7), we obtain Eq. (Fll) with Eq.(4-2). □ 

If the variations are taken in P, 5 2 E is a constant. The function -T(r) is deter- 
mined from the first variation of the metric, and thus it can be thought of as a given 
function. Applying Formula ||, we obtain 

Lemma 10 Let V(t) be the functional in the space of solutions P defined by Eq.(2-t ), 
and suppose that the spatial manifold M is closed. Then, the second variation about 
a homogeneous and isotropic solution is given by 

5 2 V(t) = a 2 a[— — + 5 2 C) -— 4-9 

where 6 2 E and 5 2 C are constants, a = a(r) is a solution of Eg. ( 3- 1 (\ ), and -T(t) is 
defined by Eq.(4-2). 

The integration constant 5 2 C appearing above is naturally identified with the 
variation of 5C appearing in Lemma ||, since their coefficients are the same. The 
constant 5 2 C is therefore the second variation of the big bang constant C. Our 
conclusion in this section is the following (cf. Theorem [?]): 

Theorem 11 The second variation ofV(r) along a smooth path in Av~ 1 o(ao, To c ) 
from an FLRW solution o(ao,TQ C ) is given by 

ev<?) = f ^dr, (440) 

2 J aa z 
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where a = a(r) is a solution of Eq.(S-lC), and -T(t) is defined by Eq.(4-2). 

§5. Comment on the boundary effect 

In the previous two sections we have mainly considered the cases of spatial 
manifolds being closed, i.e. dM = 0. Such a case is clearest from a theoretical point 
of view, since we do not need information on boundaries. However, we may also 
be interested in a compact spatial manifold with boundaries, for example, a spatial 
manifold inside the horizon. If we restrict to the case in which the boundaries are 
comoving so that the particle number E is conserved, we can apply a procedure 
similar to the ones used in the previous sections. See Lemma [l|. Now 5B does not 
vanish, but it can be shown that it is constant. In fact, the linearized momentum 
constraint equation proves 5B = 0. Thus we can again integrate Eq.(3jl) using 
Formula ||, and we obtain 



6V(t) 



a a 



SE 
~2 



dr 



aa z 



5B_ 
~2~ 



dr 



+ 5C 



(5-1) 



instead of Eq. (3-13), where 5E, 5B, and 5C are constants. 
The explicit forms of the term for the boundary effect 



nr , . 6B 2 . fdT 



(5-2) 



are the following, 
(i) For k = 1, 



SVb(t) = — —Oq(1 — cos??) (4(1 — cos??) + sinr?(sinr/ — 3r?)) . 



(5-3) 



(ii) For k = — 1, 
SV B (r) = -§ 
(hi) For k = 0, 



a^coshr? — 1) (4(coshr? — 1) — sinhr?(3r? — sinhr/)) . (5-4) 



SV b (t) 



5B_9_ 

ITIo 



ao{T - T 0c 



(5-5) 



The mode corresponding to 5Vb(t) is not obtained from the space of FLRW 
solutions. 



§6. Summary 

We have succeeded in obtaining the variations 5V(t) of the spatial volume about 
an FLRW solution in the space of solutions for a dust system (Lemma ||). It is 
noteworthy that this has been done without solving linearized Einstein equations, 
and the result contains all the cases of curvature index k. Also, we remark that we 
did not split apart 7^ = Shij into tensor, vector, and scalar parts, as is done in usual 
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cosmological perturbation theory (e.g. Ref. |13|)). Hence our result does not depend 
on such a part. 

It is, however, more surprising that only two modes (patterns of time-dependence) 
appear in 5V (r), in spite of the fact that the space of solutions is infinite dimensional. 
This implies that almost all possible modes vanish. The origin of the two modes is 
that the FLRW solution is a two-parameter solution. In fact, explicit calculations 
(see Eq.( |3T7| )) reveal that the two modes caused by variations with respect to the 
FLRW parameters (clq, tq c ) coincide with those found in the (full) space of solutions. 
Clearly, these modes are not of essential significance in themselves. Instead, their 
existence tells us that we should divide the space of solutions in such a way that in 
each piece there exists only one FLRW solution. As we have seen, this can be done 
by considering the foliation defined by the level sets E = const and C = const. If the 
variations are taken in such a set (leaf) (i.e., if 5E = and 5C = 0), the variation 
of volume SV(t) completely vanishes (Theorem ^ and ^). Each FLRW solution is 
therefore a critical point in the leaf that contains only one FLRW solution. Every 
inhomogeneous solution in such a leaf is therefore naturally related to the unique 
FLRW solution (cf. Eq.( f3-2CQ ), and this correspondence has defined our averaging 
(up to first order). 

The second variations do not vanish in general (See Eq.( 4T0| )). We cannot 



therefore define an averaging at second order in the space of solutions. One may 
need to generalize the space of solutions to define it. This remains as a future work. 
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